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Abstract 

The  existence  of  global  weak  solutions  to  an  initial  -  boundary  value 
problem  for  a  nonlinear  hyperbolic  system  which  models  fluid  flow  in  a  pipe 
is  proven.  The  effect  of  friction  is  modeled  by  adding  a  quadratic  zero 
order  term  to  the  system  of  conservation  laws  for  compressible,  frictionless 
flow.  A  priori  bounds  are  obtained  by  means  of  a  nonincreasing  functional 
that  is  compatible  with  the  friction  effects.  The  boundary  values  for  this 
problem  cannot  be  imposed  weakly,  so  new  results  on  the  regularity  of  the 
solution  at  the  boundary  are  given. 


TV 


1.      Introduction 

Fluid   flow   in   pipelines    is   usuolly  modeled  by   the  quasilinear  hyperbolic 
system 

Pt  +   G^  "^   0    ,  (x,t),  (0,l)x(0,c.)    , 

(1.1) 
^t  ^   (^^/^^x  ■'   P^^'^x  "  ~    f|G|G/2Dp      , 

where  p   is  mass   density,    G   is  momentum  density,    p  =   p(p)    is   pressure, 

f  =   f(|G|)    is   the   Moody   friction   factor,   and   D  is   pipe   diameter   [12].      In   this 

paper,  we  construct   global    weak   solutions   to    (1.1)   satisfying  given   initial 

conditions 

p(x,0)    -   Pq(x),  G(x,0)    -   Gq(x)    ,  X  r.  [0,1]     ,  (1.2) 

and  givei!  boundary  conditions 

p(0,t)  =  P^(t)  ,        t  r  (o,co)  ,  (1.3a) 

G(l,t)  =   0  .  (1.3b) 

This  poses  the  classical  "waterhammer  problem"  since  the  waterhammer 
phenomena  in  hydraulics  can  be  created  by  a  sudden  valve  closure  downstream 
(modeled  by  tfie  boundary  condition  G  ::  0)  or  by  a  rapid  change  in  the  pressure 
upstream  (modeled  by  a  discontinuity  in  p^).   These  events  create  pressure  waves 
which  are  reflected  at  the  boundaries. 

The  term  -f|G|G/2Dp  accounts  for  the  momentum  loss  due  to  viscous  friction 
between  the  fluid  and  t^le  pipe  wall.   Since  the  flow  changes  from  laminar  to 
turbulent  at  a^  flow  rate  near  G  =  2000i'/D  (where  i,  is  the  dynamic  viscosity), 
the  properties  of  f  also  change  at  G  =  G  .   In  the  laminar  region 


f(|Gl)  -  64p/1g|D  ,         IgI  <  G^  ,  (1.4) 

but  the  friction  factor  is  determined  experimentally  for  turbulent  flow 
(|g|  >  G  )  and  depends  on  the  pipe  roughness  (which  we  assume  to  be  constant 
in  space  and  time)  as  well  ason.the  flow  rate.   In  particular,  it  can  be 
observed  from  experimental  data  that  there  exists  a  constant  f,  >  0  such  that 

tun         f(|Gl)  =  f^   .  (1.5) 

I  G|->-  «> 

Thus,  the  friction  term  f|GlG/2Dp  is  nearly  quadratic  in  G  for  turbulent  flow. 
(See  [12]  for  a  discussion  of  the  theoretical  and  experimental  basis  for  the 
Moody  friction  factor).  Our  analysis  will  assume  only  the  following  properties 
for  KG)  =  f|GlG/2D  : 

H(0)  -  0  ,  (1.6) 

H,>H>0,  ^'-'^ 

H  is   locally  Lipschitz   continuous.  (1-8) 

Property   (1.6)   states  that  there  should  be  no  friction  v;hen  there   is  no  flow. 
Property   (1.7)   states  that  the   relative  change   in   the   friction   (assuming  that  p 
is   fixed)   is   greater  than  the   relative  change   in  the  flow  rate.     This   is  obviously 
valid  in  the   laminar  regime   (H(G)   =  64G/viD)   and  in   the  completely  turbulent 
regime   (H(G)   =   f,  !g|G).      Our  study  of  the  Moody   diagram  [12,   p.    297]  has   led 
us  to  assume   its   validity  in  general.      We  note  that   (1.7)    is  equivalent   to   the 
condition   (fG)p  lO   .    Property   (1.8)    is  justified  for  all    flow  rates,   G,   except 
possibly  at  the  transition  flow  rate    |g|    =   G^     (See   [7]  where   f  is   allowed   to  be 
multi-valued  at    \G\    =   G   ). 


We  also  assume   that   the  sound  speed,   c   >  0,    is   constant,    i.e.. 


P   (p)   =   c      . 


(1.9) 


This  is  valid  for  an  ideal  gas  v;hich  is  maintained  at  a  constant  temperature  by 
heat  exchange  between  the  gas,  the  pipe  wall,  and  the  surrounding  environment. 
For  many  physical  problems  property  (1.9)  is  also  a  good  approximation  for 
modeling  the  flow  of  liquids.  We  believe  that  the  perturbation  techniques 

used  in  [13]  can  be  utilized  to  obtain  results  for  more  general  equations 
of  state,  but  it  is  our  view  that  the  additional  technical  complications 
would  obscure  the  ideas  presented  here. 

In  [7],  Luskin  has  shown  for  the  initial -val ue  problem  (1.1)-(1.2)  that 
a  unique,  global  smooth  solution  exists  if  the  initial  data  are  in  an  appropriate 
invariant  region  and  if  the  first  derivatives  of  the  initial  data  are  sufficiently 
small.  However,  if  the  first  derivatives  of  the  initial  data  are   too  large, 
then  discontinuities  can  be  shown  to  occur  even  when  the  data  is  smooth  (This 
can  be  done  using  a  variant  of  Lax's  ideas  for  the  frictionless  case  [4]). 
To  allow  for  more  general  data  here,  we  need  to  consider  weak  solutions  of 
(1.1).  We  call  p,GcL"(q)  a  weak  solution  of  (1.1)  if 


j  [P^^   +  G(^^]dxdt  =  0  , 


n 


[Gi>^  +    (G_     +   p(p))^^  -    (f|G|G/2Dp)<^]dxdt  =   0      , 


for  all  (^  c  C^   (n),   where  n  =   (0,l)x(0,<»). 

The  principal  result  of  this  paper  is  the  following  theorem: 

Theorem  1 .  Assume  that  properties  (1.6)-(1.9)  hold  and  that 


(1.10) 


Var  In  p^  +  Var  £n  p„  +  Var  ^    <    In  ^-^     -    .96        .      (1.11) 

t>0         ^       xe[0,l]  °       xc[0,l]   ^^0  "^ 

oo 

Then  there  exists  a  weak  solution  p,  G  e  L  {n)   to  (1.1).   The  initial  values 
are  satisfied  in  the  sense  that 

p(-,t),  G(-,t)e  Lip([0,"),  J(O.l))  (1.12) 


and 


tim     p(-,t)  =  p^   ,     JUm     G(-,t)  =  G.  in  L^O.l) 
t-^o        °   t^O        ° 


The  boundary  values  are  satisfied  in  the  sense  that  for  any  T  >  0 

p(x,.),  G(x.-)e  Lip  ([0,1],  L^O,T))  ,  (1.13) 

and 

-ton  p(x,-)  =  pg  ,  Lim     G(x,-)  =  0  in  L  (0,T)  . 
X  ^0  X  ->1 

(Here,  e.g.,  p(-,t)£  Lip  ([0,T],  L  (0,1))  means  that  there  exists  a  constant,  C, 
such  that 

|p(-,tn)  -  p(-,t2)i  1     <  C  I  t.  -  tp  I  , 
'  L  (0,1) 

for  all    t, ,   t2  e[0,co)).     Without  loss  of  generality  we  assume  that  Pq(0)   =  pd(0) 
and  that  6^(1)   =  O  by  redefining  Pq(0)   and  6^(1)    if  necessary.      In   this 
way  the  incompatibility  of  initial   and  boundary  data  is  accounted  for  in 
the  left  hand  side  of  (1.11)   by  allowing 

tim     P    (x)  ^  p   (O)   and  Urn     G  (x)   ^   Q  . 
x-^O     °  °  x-^1 

The  only  purpose  of  (1.11)  is  to  guarantee  a  priori  that  the  flow  remains 
subsonic;  i.e., 

Ivl  <  c  for  (x,t)E  n   ,  (1.14) 


where   v  =   G/p   is   the   velocity  of  the   flow.      Note   that    (1.11)    is    independent 

of  units,   and  is   a  satisfactory  bound  for  subsonic   flow  since    |v|  <c  is   sustained 

by  a  bound  of  as   much   as    .96c  on   the   variation    in   v  and  c  £iip   of   the   data.  The   a  prio 
bound   (1.14)    is   required   to   guarantee   that   the  boundary   conditions    (1.3)    can 
be   imposed.      In   general,    boundary  value  problems   for   (1.1)    in  which  cither  the 
density  or  the   flow  rate   arc  assigned  at  each   boundary   can  be   solved   uniquely  only 
when   the  characteristic  speeds   A-j,    A2  satisfy   a-j    <  0,    a^   >  0   .      Our  problein  is 
posed  in  Eulerian   coordinates  where  the   characteristic  speeds   are   A-i    =   v-c, 
Ao  =   v+c   ;  so   (1.14)    is   required  for  A,    <  0    ,    A^   >  0  .      Earlier  work  on   the 
construction  of  solutions   to   initial -boundary  value  problems   has   been   done 
by  Nishida  and  Smoller   [11]   and   Liu   [S]   for   the   "piston   problem",    but  a  priori 
bounds   similar  to    (1.14)   were  not   required   there.      This   is   because   the   piston 
problem  is   posed  in  Lagrangian   coordinates  where   the  boundaries  move  with   the 
fluid.      Thus,    for  the  piston   problem,    ^-j    =   -c,    ^2  ^   ^'   ^0   ^i    <  0,    A^   >  0   is 
already  guaranteed  a   priori. 

Note  also  that  boundary  condition   (1.3b)    is   a   "natural"   boundary  condition 
and  could  have  been   imposed  weakly  by  requiring  that 


n 


[p4,^  +  G6^]dxdt  -  0 


for  all    (f- E  C"  ((0,1  ]  X   (0,^o)).      However,    the  boundary  condition    (1.3a)    is   not  a 
natural    boundary  condition   and   it   is  necessary   for  us  to  give  new   results  on   the 
regularity  of  the   solution   at   the   boundary   in   order  to  make  sense  of  boundary 
condition    (1.3a).      This   problem,   as  well,    did   not   arise   in   [5]  or   [11]   since 
the   boundary  conditions   for   the   piston   problem  are   "natural"   boundary  conditions 
and  could  be   imposed  weakly. 

Our  method  is  to  obtain  a  solution  to  (1.1)-(1.3)  as  the  limit  of  approximate 
solutions  which  are  constructed  by  a  fractional  step  procedure.  In  the  first  part 
of  Oc'.ch   step  we   use   Climm's   method   [1]   to  approximate   the   solution   of  tno   system 


of  conservation  laws  for  frictionless  flow.  The  second  part  of  each  step  accounts 
for  the  effect  of  friction  on  the  flow^  and  involves  solving  an  O.D.E.  that  is 
determined  by  the  zero  order  term.  Liu  [6]  and  Ying  and  Wang  [14]  have  also 
proven  the  existence  of  global  solutions  for  some  systems  of  conservation  laws 
with  zero  order  terms;  and  the  latter  used  a  fractional  step  method  similar  to 
ours.  However,  their  analyses  took  account  of  only  the  magnitude  and  not  the 
orientation  of  the  vector  field  given  by  the  zero  order  terms.  These  methods  are 
inadequate  for  our  purposes  because  the  physical  friction  term  fJG|G/2Dp  is 
quadratic  in  G  at  infinity;  and  solutions  to  systems  of  conservation  laws  with 
quadratic  zero  order  terms  will  "blow  up"  in  finite  time  if  the  associated  vector 
field  is  allowed  to  have  an  arbitrary  orientation.  Moreover,  the  methods  in  [5] 
"Bnd  [14]  will  not  imply  the  a  priori  bound  (1.14)  unless  the  orientation  of  the 
vector  field  is  considered.  Thus,  because  our  friction  term  is  quadratic  in  6 
and  because  our  boundary  value  problem  requires  an  a  priori  bound  for  the  solution, 
it  is  crucial  that  we  found  a  nonlinear  functional  which  is  equivalent  to  the 
variation  norm  and  which  is  nonincreasing  on  both  of  the  fractional  steps. 
Although  the  functional  introduced  by  Nishida  [9]  is  nonincreasing  for  the  system 
of  conservation  laws,  it  is  inadequate  for  our  purposes  since  it  can  increase  on 
the  friction  step.  However,  we  show  that  if  the  zero  order  term  satisfies  certain 
monotonicity  conditions,  then  the  functional  given  by  Liu  [5]  is  nonincreasing 
for  both  fractfonal  steps.  These  monotonicity  conditions  are  satisfied  by  the 
physical  friction  term  in  (1.1)  when  the  flow  is  subsonic.  The  fractional  step 
procedure  that  we  have  analyzed  has  recently  been  developed  by  Marchesin  and 
Paes-Leme  [8]  to  obtain  numerical  results  for  system  (1.1). 

In  section  2,  we  discuss  the  solution  of  Riemann  problems.  The  fractional 
step  scheme  is  defined  in  section  3  and  the  basic  stability  result  is  given  there. 
In  section  4,  the  structure  of  Riemann  problem  solutions  is  further  analyzed  to 
obtain  bounds  on  the  "nonlinear  interaction";  these  bounds  are  used  in  section  5 


to  prove   regularity   results   for  the   approximate   solutions.      Section   6  contains 
the  analysis   of  the  convergence  of  tfie  approximate   solutions. 


2.      Solution  of  the  Riemann    Problem 

The  Riemann  problem  is   the   initial    value  problem  for  data  which 
is  constant  to   the  left  and  right  of  x  =  0.     We  study   the   Riemann 
problem  for  the  nonlinear  hyperbolic  system 

u^  +   F(u)^  =   0,  (x,t)eRx  R"^        ,  (2.1) 

where  u  =   (p,G)^'"   .   F(u)   =   (G.   G^/p  +  pCp))^"",   and  p'{p)   =   c^   .      The 
eigenvalues  of  dF  are 

A^(u)   =  V  -  c   .  \^{u)  =  V  +  c  (2.2) 

with  corresponding  right  eigenvectors 

R^Cu)   =   (1.   V  -   c)*"^  ,  R2(u)   =   (1,   V  +-  c)^^     . 

The  main  existence  result  is  that  for  initial   data 

EU|      if     X  <  0     , 
Uf^     if     X  >  0 

(we  always   assume  p.  ,   pp  >  0)   there  exists  a   unique  solution   u(x,t)   =  u(x/t) 
such  that  u(x/t)   consists  of  constant  states   separated  by  "shock  wave"   and 
"rarefaction  wave"   solutions   [9]. 

We   first  discuss   the   rarefaction  wave  solutions.      We  note  that 
a  smooth  solution  \j{K),K=  x/t,  must  satisfy 

[dF  -  5l]u(0   =   0       ■ 

Hence,  a  smooth  solution  o{e,)  must  satisfy  u(£;)e  span   {r^{u{e,)))   and 
5  =   X   (u(c))fo'')i  =   1    or  £   =   2.      An  £-rarefaction  wave   is   a   continuous 
solution,   u(x,t),  whose   values   lie  on   an   integral   curve  of  the  eigenvector 


V     .     The  functions 


c    :=    V  +    c  lirp      ,      „   ,    y   -_  ^  ^M^    , 
2  ^  -  2     '      ' 

are   Rieniann   invariants;    i.e., 

v^   s-R^    =  0   ,  \r-^~2  =  0     , 

and  s[resp.  r]  is  constant  on  an  integral  curve  of  R-.    [resp.  Rp].   Thus, 

the  ^-rarefaction  curves  can  be  defined  by 

RjC^l)  =  fup  I  r(uj^)  >  r(U|_),  s(uj^)  =  s(U|^))   ,       (2.5a) 

R2(U|^)  =  {uj^  |r(up)  -  r{u^),    s(up,)  >  s(U|^)}  (2.5b) 

=  (Up  |V|^  -  Vj^  =  -czforz  =  -^ipn  -  -^^'Pi  i  0) 

It   is   important   to  note  that 


^:    (P,G)   -   (r,s) 


is  a   1-1    regular  map  of  R     x  R  onto  R  x  R 


A  1-shock  viave   [resp.    2-shock  v;ave]  of  speed  o   is   a  weak   solution 

u(x,t)   -     {    ^  (2.6) 


u        i  f     x/t   >  o  , 

K 


which  satisfies  the  Lax  entropy  condition  [3] 


[resp. 


>^1  (u^)  >  o  >  A^  (uf^) 


^^(ul)  ^  o  >  >'2(uj^)]. 


(2.7a) 


(2.7b) 


Since   u   is   a  weak   solution,    it  must  satisfy   the  Rankine-Hugoniot  jump 
condi  t  ion 


o[u^_  -   Uj^]  -  F(u^^)   -   F(uj^)      .  (2.8) 

By  eliminating  a  in   2.8)   and  applying  the  Lax  entropy  condition  we  obtain 
the  following  £-shock  wave  curves 

V^l)   =    (Uf^:   v^-Vj^  =   c(e~'/2  -  e  '/^)  (2.9a) 

for    2  =  liip.    -  Inp^  ^  0}    , 

S^Cu^)   =    {u^:   v^-Vj^  =   c(e-^/2-e'/2) 

for    z  =  £j:p[^-£np|_  <  0}    .  (2.9b) 

Substituting  in  (2.8)   gives 

o  =  Vl  -  c  e"^^^  =  v,^  -  c  e^^^,  z  =  -^p^  -  -^"Pj^  1  0    ,        (2.10) 

for  a  1-shock  and 

o  =   V.    +  c  e"^/^  =   Vj^  +  c  e^^^    ,  z  =  Inp^  -   Inp^   <_  0     , 

for  a  2-shock. 


It  now  easily  follows  from  (2.5)  and  (2.9)  that  z  ^lup.-   Inp 


R 


[resp.    z  =  ^ipp  -  £np.  ]   is  a   regular  parametrization  of  the  C     curve  T-j(u^)   = 
R^(u,_)US^(Ul)    [resp.    T^    (u^)   =  R2(ul)US2(ul)].     We   call    z   the 
"signed  strength"  of  a  given  wave   (so  that  the  signed  strength  of  a   rare- 
faction wave  is  positive  and  the  signed  strength  of  a  shock  wave  is  negative), 
and  we  call    |z|    the  strength  of  a  wave.     The  existence   theorem  for  Riemann 
problems   follows  directly  from  the   fact  that  given   any  two  states   u^  and 
Up,   there  exists   a  unique  state  u^  such   that  u^^  e  T^(U|^)   and  Up^eT2(u^) 
[9]-,   i.e.,    the  Riemann  problem  for  (2.1)    can  always   be   uniquely  solved  by 
a  1-wave  that  connects   u,    to  Uj^.  and  a  2-wave   that  connects   u^  to  u^  . 

It  is   useful    to  view  the  rarefaction   and  shock   curves   in   the  rs-plane  . 
Since  r[resp.    s]   is   a   Riemann   invariant,    the  1-rarefaction   curves   [resp. 
2-rarefaction   curves]  are  parallel    to   the   r-axis   [resp.    s-axis]. 

10 


Also,  the  1 -shock  curve  (?.9a)  is  easily  shown  to  be 

^L  '  ^R  "  f  ^^'^^^  -  e  ^''^  -  z),  2  =  £);n^  -  Inv^   <  0  ,     (2.11a) 

^L        ^R       2    ^  ^ 

and  the  2-shock  curve   (2.9b)    is   easily  shown   to  be 

Tl   -    r^^  =   2'  (^"^^^   -   e   ^^^  +  ?),      •      z  =  £npp  -  ^np^  i  0   >  (2.1  lb) 

c    ,    -z/2  z/2  > 

s^-   s^=  j{e  -  e  -   z). 

We  can  obtain    from   (2.11)    the   following   theorem  concerning   the   geometry 

of  shock   curves   in   the  rs-plane. 

Theorem  2. 1    [9]     All      -^-shock   curves   of  equal    strx'iigth  are   translates 

of  one  another.      The   2-shock   curve,   S2(u,  ),    is   the   reflection  of  the   1-shock 

curve,   Si(u.  ),   with   respect   to    the   line  ^V^P   =   Inp      .      Morever,    if  Uj,  =   Un(z), 

z   <_  0,   denotes   the   parameterization  of  ^1(^1  )    given   above,    then 

and      -r      5(ud(z))/   -r-     r(ur,(z))    is   monotone  decreasing  from  1-   at  -co     to 
dz  K  az  K 

0+  at  0   . 

Finally,  we  shall  need  to  construct  the  solutions  to  certain  initial- 
boundary  value  problems.  When  ttie  boundary  is  x  -  0,  we  consider  the 
problem 

u^  +  F(u)^  =  0  ,  (x,t)t:  R"^  X  pJ      , 

u(x,0)  =  Uq(x)  =  Up^  ,     X  i-  R  (2.12) 


p(C,t)  =  p,    ,  tcR' 
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It  can  be  checked  that  there  exists  G,    such  that  u     e  T„(u, ).     But  the 

L  K  c      L 

2-wave  connecting  u.  to  Up,  will  take  the  value  u.  at  x  =  0  only  if  it 

has  positive  speed.  If  u^  e  Rp(u.  ),  then  it  is  necessary  that 

Xp(u.  )  =   V  +  c  >  0  to  f^uarantee  that  the  2-wave  has  positive  speed.  If 

Ud  e  S  (u,  ),  then  the  corresponding  2-wave  has  positive  speed  if  v^  >  -c 

since 

o  >   X2{Uf^)   =      v^  +  c  >  0   . 
When  the  boundary  is  x  =  1,  we  consider  the  problem 

u^  +   F(u)^  =  0   ,  (x,t)   e  (  —  ,1)   x  R"^      , 

u(x,0)   -   \    ,       ■  X  e    (-00,1),  (2.13) 

G(0,t)   =  0   ,  t  £  R"^       . 

tr 
In  this   case   there  exists    Pn  such   that  U[^  =    (pp,0)       e  Ti(u^),   but  the 

1-wave  connecting  u.    to  Un  will    take   the  value  Un  at  x  =   1   only  if  it 

has  negative  speed.      This   is   true  if      v.    <  c.     Thus,   the   initial -boundary 

value  problems    (2.12)   and   (2.13)   can  be  solved  by  simple  waves   if  all 

the  velocities  occuring  in  the  solution   are  subsonic. 
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3.      Definition  and  Sta[)ility  of  the   Fractional    Step  Scheme 

In   this   section,   v;c   construct  approximate   solutions   to   (1.1)-(1.3). 
Let  h  -^    ]/N,   N   a   positive   integer,   x^    =    ih,I.    -    rx_i>   x.],    1   =    [0,1], 
and  let  k  >    0,   tj   -  jk    ,      J.   =   [t.-j  ,    t  •].      Also,   let   a  =    (a^  ,a„, . . , ) 
a.  t:(0,l),   he  a   sequence.      We   define  approximate  solutions   u,    =    ( p^ ,    CkJ 


tr 


h        ^^^h'    ''U' 


tr 


and  u.    =    (pu,    G^)        inductively.      Assuine    that   u,    and  u,    are   defined   for 

t   <  t-    .      Then   u^,  on    IxJ.,,    is   the   solution    to 
-    J  h  j+l 


Uu,    +    f(iii  )      =   0    , 
h  t         Mr  X  ' 


(x,t)  e    IxL^^ 


p\(t,0)    =   P,(t.,J, 


t   c    J. 


(3.1) 


G^(t,l)   =   0      , 


t  c    J 


J^l       ' 


Uh(x,t..)    .    u^,(x._^     ^    n.h,    t.-), 


u^(x,00    ^    u^(x._J    , 


X  c  I .  ,    i  f  j    >  0    , 


x  c  I.,    if  j    -   0 


where   F(u)   =    (G,    G^  /p     +    p(fO)^^  and  u^(x)    =    (p^(x),    G^(x))^^ 

0  V  0    ^  ^         0    ^ 

Next   define   the   functions   u(t,u)    =    (p(t,p,G),   G(t,p,G))   by 


u^  =   H(ti)    , 
u(0)    --=   u 


t  >  0   , 


(3.?) 


where  H(u)  =  (0,  -f  |  G|  G/2IV.  )^'"  -  (0,- H(  G)/p  )^^  .  Then  we  set 


u^|(x,t)  -  u(t-t.,  u^^(x,t)),      t  €    J.^-| 


(3.3) 


Note   that   (3.1)    poses   an   initial    value   Riernann   problem  at  each  mesh 
[loint   (x.,t.)   0  <   i    <   N,   and  a   boundary   Rieinann   problem  of  ty[)L'   {?.]?) 
and   (?.13)    at    i    ^   0  and   i    -■=   fj,    res[)ecti  vely .      Therefore,   we   can   use   the 
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Riemann  problem  solutions  of  section  2  to  solve  (3.1)  in  J-,,,  as  long. 

as  the  waves  in  these  solutions  do  not  intersect  in  J-j+t  >  and  as  long 

as  |v|  remains  less  than  c  in  the  boundary  problems.  Given  these  two 

conditions,  Ul,  and  u.  are  defined  for  all  time.   Choose  ^  fixed  with 
n     h  k 

^>  4c   . 

We   guarantee   the   two  conditions   above  by  showing  that   if  (1.11)   holds, 
then    |v|    <  c   for  the  approximate  solutions   u,  ,   and  u.     .      It  then   follows 
from  (2.2)  that  the  wave  speeds   are  bounded  by  2c,   and  so   the  waves   cannot 
interact  in  J-^-.    . 

We  shall    use  "I-curves"    in  our  proof  of  the  stability  of  the  scheme   [1]. 
We  call   the  points 

(x.  +  cxj.h,   tj),  0  <  i   <  N-1,     j  ^0   , 

(x.    ,   t.,,  )    ,  0  <   i    <  N    ,       j   >  0   , 

vertex  points.      An   I-curve,   J,   is  a   curve   in  xt-space  which   successively 
connects   vertex  points  on   tj  to  adjacent  vertex  points  on   t-    j^      such 
that  the  index  i    is  a  nondecreasing  function  of  x  on  J,  and  such  that  J 
is   linear  between  adjacent  vertex  points   (see  Fig.    3.1). 
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t  =  2k 


t  =   k 


t  =   0 


X  =   0 


X  ■-   1 


->   ^ 


Fig.    3.1 
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It  is   important  to  note   that  the  unique   I-curve,   J(j),   that  connects  the 
vertex  points  on  t,  to  the  vertex  points  on   t .  ,    crosses  all   the  waves 
in  the   Riemann  problem  solutions  of  u,    in   J.   ,    .     We  partially  order 
the  I-curves   by  saying  that  larger  curves   lie  toward  larger  time,  and 
we  call   Jp  an   immediate  successor  of  J     if  J     connect  the  same  vertex 
points  except  one  and  if  J„  >  J,    .     We  also   let  t(J)    denote   the  value 
of  t  where  J  intersects  the  line  x  =  0  .      For  any  I-curve,  J,  define 

L(J)   =  IIyaI      , 
J     ^ 

^  t     >   t(J)  ^     ^    '  ^     ^"'' 

J 

where  L     is  the  sum  over  all  waves  which  cross  J  and  where  y.  denotes  the 

strength  of  the  wave. 
Lemma  3.1  Suppose  that 

L(J(j))  +  Lg(J(j))  =  Vj  <  In^ 


(3.4) 


Then 


and 


sup 
'x.tJE  Ix. 


In 


P^(x,t) 


-Vj 


J+1 


sup 

(x,t)G  IXjj  +  -, 


v^(x,t) 


v. 


<.^-^  <  1 


V./2 


(3.5) 


(3.6) 


where 


p_^  ^  iA.m     Pg(t) 


(3.7) 


Proof.  The  proof  of  (3.5)  is  an  immediate  consequence  of  our  definition 
of  wave  strength.  To  prove  (3.6)  note  that  v^(l,t)  =  0  and  that 
u.(l,t)  is  connected  to  U|^(x,t),  t  e  J.^.-|  ,  by  a  curve  which  is  composed 
of  a  sequence  of  shock  and  rarefaction  curves  the  sum  of  whose  strengths 
is  bounded  by  V-  .   By  viewing  -^-waves  in  the  rs-pTane,  it  is  evident  from 
Theorem  2.1  that  |v^(x,t)|  must  be  smaller  than  the  change  in  v  across 
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t  r 
an  ^-shock  which  starts  at  u  =  (n  ,0)   and  has  strength  V.  .   For  suet 

a  shock  v;avc,  we  can  use  (?.lla)  to  calculate  that 


V. 
e  J  -  1 

^  J 


We  next  prove 

Theorem  3.1      Suppose   that   (1.11)   holds, 


Then 


V  7— 

V     <   V   H   Var     Inp^  +  Var     Inp     +    Var     ^     <  In  -^^-^ 
0  -  X   .    r,       B       „   ro  1 1  0    V.  rn  n  c  2 


t  >  0 


xc[0,l]  °    xl[0,1] 


(3.8) 


and 


L    (J(j  +  1))    +   L(J(j  +  l))<   L    (J(j))   +   L(J(j)) 


(3.9) 


for  j   =  0,    1 ,    ...    , 

The    inequality    (3.8)    follows   easily   from  Theorein  2.1,   so   that   the 
remainder  of  the   theorem   is   a   consequence  of   (3.9). 
We  shall    show  that 


L    (J2)   +   1(0^)    <  L^(J^)   +   L(J^) 


(3.10) 


whenever  Jo  is  an  immediate  successor  of  J-,  .  To  verify  (3.10),  it  suffices 
to  study  the  "interaction"  that  occurs  in  the  region,  a,  that  lies  between 
Ji  and  J2  ;  i.e.,  it  suffices  to  compare  the  strength  of  the  waves  that 
cross  J^/J^  to  the  strength  of  th(^  waves  that  cross  J^/J2  •   If  A  does  not 
contain  a  mesh  point  (x-,  t-),  then 

Lg(J2)  +  l{0^)    -   L[.,(J^)  +  L(J,)   . 
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Thus,  we  need  only  consider  the  case  when  A   is  a   "diamond"   centered  at 

(x-,   t-)   for  1   j<   i   £  N-1    or  a  "half-diamond"   centered  at   (x   ,    t.)   or 

/ 
(xm,   t-)-      We   denote   these  diamonds   by  a-^    .      We  let  y.-   denote  the  signed 

strength  of  the  £-wave   in   the  solution  of  the  Riemann   problem  which   is 

posed  at   (x-,    t.)    .      Here,   y..   denotes  both  the  name   as  well    as   the  signed 

strength  of  a  wave. 
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3+h 


G^(X._^     .    a.h,t.:0 


J       2 


"r       %^'^^    '   "j^'S"^ 


^■-1   -^  "j^' 


X.    +    a-h 
1  J 


Aij    {~^dU      =    <i> 


'j^H  ^B^^jW    =  ^ 


Uf^  :    u^(-jh,tj') 


>^o   =   0 


X   =   a  .h 
J 


A.j  n    {x=0}   /    <■ 


^j^^. 


'3-H 


%^\-]  '  "j^'^Y^    ^ 


^N   ^  "jh  ^N   ^'    ^ 


A..n(x=^i}  /  <• 
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Since  our  scheme  has   two  fractional    steps,  vie  will    need   to  define 

L(J2)  ""  ~      Iy/\I    to  be   the   sum  over  all   waves  which  cross  02^10,   plus  the 

sum  of  the  waves  which  solve  the  Riemann  problem  with  constant  states 

u[  and  Up  (See   Fig.    3.2).     We  also  define   L   (J2)   =  •-    (J2)    . 

Lemma  3.2     (Liu,    [5])      If  J2   is   an   immediate  successor  of  J-,,   then 

Lg(32)   +  L(J2)   <  Lg(J^)  +  L(J^)      . 

The  next  lemma  says  that  the  functional   L     +  L  is  also  nonincreasing 

B 

for  the  friction  step. 

Lemma  3.  3  If  Jo  is  an  immediate  successor  of  J-,  ,  then 

Lg(J2)  +  HJ^)  1  ^^C^2^   +  L(J2)  •  (3.11) 

The  proof  of  Lemma  3.3  is  given  in  Appendix  A.  Now  (3.10)  follows 
from  Lemma  3.2  and  Lemma  3.3.   It  follows  from  (3.10)  that  if  J^  >  J-,, 
then 

Lg(J2)+  L(J2)  1  Lg(J^)  +  L(J^)   .  (3.12) 

The  stability  result  (3.9)  is  a  special  case  of  (3.12)  since  J(j+1)  >  0(j). 
Finally,  we  note  that  the  stability  results  in  this  section  did  not  require 
any  continuity  properties  for  H(G).   In  Appendix  A,  it  is  shown  that  the 
properties  (1.6)  and  (1.7)  are  sufficient  when  H   is  extended  to  be  a  maxi- 
mal nwnotone  function.  However,  the  proof  of  convergence  for  our  scheme 
will  require  the  further  assumption  that  H  be  locally  Lipschitz  continuous. 
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4.      Bounds    for  the   lotal    ^'ot)^  incari  ty 

In   this   section,   vie   shall    give   bounds   for  the   total    nonl  ineari  ty. 

These   bounds   are   required   for  the   regularity   results    in    the  next   section 

3+  'b 
It   follov/s    from  Lemma   3.1    and  Theorem  3.1    that   if  V   <  In    -^--   ,    then 


e^-l 


and 


\ln 


V/2     ^ 


<    V 


Hence,  it  follows  that  there  exists  a  compact  set,  KC(0,o:.)x 
of  h,  k,  and  a,  such  that 


(4.2) 
,  independent 


Uj^(x,t)£:  K   for  (x,t)c(0,l)x(0,")  .  (4.3) 

In  the  following,  let  C  denote  a  generic  constant  (dependent  on  c,K,  and  V) 
which  is  independent  of  h,  k,  and  a  . 


Recall  that  y^.  ■  is  the  signed  strength  of  tfie  i'-wave  in  the  solutior;  of 

f 
n  problem  for  u,  at  (x-,  t-)  and  let  ,•  :.  denote  the  sum  of  the  signed 


Ri  eman 


strengths   of  £-waves  v;hich  enter  a—    .      We   define   the   rionl  i  neari  ty  of  the 
interaction   at   (x-,    t.),   N.-,    to  be 

N.^.   =    /.   -    y:^-.    ,      0   <   i    <  N,   j    >  0    , 
ij         'ij         'ij 


OJ 


^Nj 


Y 


I  =    1 

) 

2 
oj 

2       ^       ■> 

"^   OJ         ^  OJ 

I    =   ? 

» 

1 
'Nj 

5 

1            2 

1 
2 
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"fj 

.y4.  , 

IJ    ' 

%^J 

4-V 

V    OJ 

OJ            0 

Similarly,    if  j'^       is  the  signed  strength  of  the   £-wave  that  crosses 
J^/J-,    ,   we   denote   the  nonlinearity,   N..      ,   that  occurs   in   the   conservation 
law  step  by 


0  <  i    <  N     . 

£  =  1      , 
^  -  2      , 


y  Mi    ~  Y  M-j    ~  Y  M-;     >  I 


£  =    2 


Finally,  we  set 

U         '   ij  ij 

N..   =  I  n!.  1     +    1  N^.  1 
IJ       '     ij  '  '     ij  ' 

We  shall:  prove  the  following  theorem. 

Theorem  4.1   For  any  T  >  0,  there  exists  C  =  C(T)  such  that 


ll         -1         -2 

Nj  "  ^  Hi  "  "^  rd  ' 


N{T)    =       I         N..     ^C(T) 
t.    <  T     ^^ 


(4-4) 


where  the  sum  is  over  (i,j)   such  that  t.  ^  T. 

Two  waves    Y    and    Y    that  cross  J  are  said  to  approach   if  they  are  both 
of  the  same  family  and  if  Y/\  °'"  Yg  is  a  shock  wave. 
Define 

Q(J)  =     I  IyaMybI  (4.5) 

App(J) 

where  the  sum  is  over  all  pairs  of  approaching  waves  that  cross  J.  Also,  let 


Dij  =    y       \y;^\\\\  (4-6) 

'^  App  (J^/J^)  ^   ^ 
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be. the   sum  of   the  products  of  approachinn  waves   that  enter  A..      .      for 
a   simple  wave    y  ,   we   define   the   sftock  wave   strength   and  rarefaction 
wave   strength  ofy    by 


S(y) 


R(r) 


Then  we  decompose  D..  into 


1     ? 

D.  ■  -  D..  +  0  , 
1J    IJ     iJ 


whore 


Dij  ^  I      S(Y,)R(Yp) 
App(J^\J^) 


and 


Dfj  =  I     S(,,)S(.^) 


^2 

i 

App(J^\Jp) 
We  shall  need  the  following  lemmas  in  the  proof  of  Theorem  '1 . 1 


Lemma  4.1  There  exists  a  constant,  C,  such  that 

^j  -  ^^ij  '         0  <  i  <  N  ,  (4.6) 

\^   1  |aL|   ,  (4.8) 

where  AL   =   KJ^)    -   L(J^)    and  aI  ^   -   L^,    (J^)    -   Lp    (J^)      . 
Proof.      Note   that   if  a   single   1-wave   and  a   single   ?-wave   enter  A..    ,    then 
N^  .   =   0   .      Thus,    (4.6)    follows   from  the   results   in    [10,    13].      The  estim-iles 
(4.7)   -    (4.8)    are  easily   verified  by   checking  all    possible   boundary   inttn-- 
actions   for   the   conservation   law  step   [S].  El 
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Now 


^^^^ij  ~~    (^  '  ^ij)  ^"^  ^j=^^]j  '  ^?j^  ^^^'"^ 


T,j  =  (0,  Y^)  andT,.  =  (yl-^   0)). 
Lemma  4.2  There  exists  a  constant,  C,  such  that 

'^ij  -  ^ijl^Ck.lY.jl   . 

Proof.   First,  let  0  <  i  <  N  .   It  follows  from  the  results  in  Section  2 

2  2    ? 

that  there  exists  a  regular,  C  ,  diffeomorphism,  6:1  -*  R  ,  such  that 

for  constant  states  u,  and  Un  we  have 

y=   eCvCu^)  -  'J'(u^))  (4.9) 

1   2 
where  Y-  (y  ,y  )  is  the  vector  of  signed  strengths  of  waves  in  the  Riemann 

problem  [u,  ,  u^l  and  where  '^((p,G))  =  (r,s).  Now  since  .Y(a)  =  '^' (fj)H(Jj) , 

L    K  t 

we  obtain  ^ 

|[4'(^(t,u[))  -  4'(u(t,UR))]  -  [>^(u[)  -  y{ur)]|  (4.10) 

<  Ct|4'(u[)  -  ^(U^)! 

Thus,  since  u^  =  u(k,u[)  and  u^  =  u(k,  u"^)  we  have  that 

\lHul)   -   y(u^)]  -  [y(u[)  -  Hu'^m  (4.11) 

<  Ck|4'(u[)  -  4'(Up)| 

Hence,  we  obtain  from  (4.9)  and  (4.11)  the  estimate 

I  ^ij  -  ^ijl  1  CI[t(u^)  -  ^(u^)]  -  [4'(u[)  -  y(u"^)]|  (4.12) 

<  Ck|Y(u[)  -  T(u^)| 

lCkI^..|    . 

This  is  the  result  for  0  <  i  <  N.  An  easy  check  establishes  the  result 
for  i  -  0,  N.  H 
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(4.13} 


We  now  define  two  nonlinear  funcLionals,  F  (J)  and  F  (J),  such 
that  r,,,(J)  decreases  by  C-jN-  •  between  successors  J-,  and  J^^  whenever  a- ^ 
does  not  intersect  the  boundary  that  lies  opposite  x  =  m.   First,  define 

i.e.,  L  (J)  is  the  sum  of  the  sti-engttis  of  all  £ -waves  that  cr-oss  J  [->]. 

Next,  define 

F^(J)  =   L(J)  T  2Lg(J)  +  r'^L-,(J)  +  ^2  Q(J)  , 

F^(J)  =  L(J)  +  M^L^CJ)  +  M2Q(J)   , 

where  M,  and  M„  are  positive  constants  which  will  be  chosen  later.  We 
note  that  by  Theorem  3.1 

F^(J)  1  C  .  (4.1b) 


(4.14) 


We  adopt  the  notation 

aL^  =  1-£(J2)  -  Lj(J^)  ,        £  =  B,  1,  2,  or  absent, 
aQ  =  Q(J2)  -  Q(Ji)  , 

as'  -   S(r^)  -  Si^.)    ,  I  =   ^..    2      , 

Ar'  =  R(y'j)  ~  R(?^ij)   .       £=1,2, 

where  S(y^-)      [resp.    R(y-^'-)    ]   denotes  the.  sum  of  £-shock  wave   strengths 

[resp.  J^ -rarefaction  wave   strengths]  which   enter  A- •    .      The   following 

three   lemmas   give  bounds   for   the   above  quantities. 

Lemma  4. 3.      There   is   a   positive   constant,  £=5(K)  >    0,   such  that   for  0   <   i    <  N, 

AL  1  -6AS^    ,  £  =   1,2.  (4.16) 

Proof.   See  [13].  C3 

Lemma  4.4.   I  f  0  '-  i  <  N ,  then 
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Proof.  See  Appendix  B  .  B 


£  =  1.2 


C4.17) 


Lemma  4.5   The  following  estimates  hold: 


I^V^V^I'^'l  i  l"-ll  *    i"-B 


Proof.      These  estimates   are  easily  checked  by  using  Theorem  2.1    [5]. 


(4.18) 


The  following  is   the  crucial    lemma  in  our  proof  of  Theorem  4.1 
Lemma  4.6 
Set 


M. 


^. 


^2  1  "• 


(4.19) 


and 


M2  1     6/2(V+l)      , 


for  c  satisfying  V  <  ^  <  1    .     Then  there  exists  a  positive  constant,   C-j , 
such   that 


'J'2^  -  '.^'^^  ^   -S^J 


(4.20) 


whenever  A.  ■  does  not  intersect  the  boundary  that  lies  opposite  x  =  m 
Proof.  We  do  case  m  =  0  (the  argument  for  m  =^  1  is  similar).   First, 
assume  that  i  >  0  (A.-Dan  =0). 
In  this  case,  we  have 

ALg  =  0  , 
AL^  -  AR^  +  AS^   , 
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aQ   =        I        IyIIas""   +   S{yl)iM'^ 

1J  TJ 


(1.21) 


;ince    IaR^I  <  of .  /2  and       I  i  Yn  1    <  V  <  1    . 


j^nJ2 


aQ 


<     I  U]^|   AS^    .     I  |,2|    ,s2  -    d]      -   Df 


J^OJ^ 


j^nj^ 


/2 


(4.22) 


Thus, 


1  ,,    J 


AFp   <  aL  +   {M^    +  M^     I  Iy^IIaS 

+    {Mo  I      Iy1|}aS^-    M,    dL   +   M,aR^    -   M„   D^./2 

^j^nj^  ^   ij      1         2   ij 

Now  since  Vi^   £  6/2(V+l),  Lemma  4.3  implies  that 

But  M,/M„  =  ^<1,  so  Lemma  4.4  implies  that 

^''o  -  ■  ^2  f^ij  -  [M2(1-C)/2]D^ 


1  -  [M2(l-0/2]D.j 

The  result  (4.20)  now  follows  immediately  from  Lemma  4.1. 

Next,  we  consider  the  case  i  =  0  (when  A.  .  H  3S^  f   j5).  Then  by 
Lemma  4.5 


(4.23) 


AL  +  2AL„  <  Ll^  <  0   , 
AL^  <_  0  , 


(4.24) 
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aQ  1    ^        IyJ   |ay  I 


—  J-]  oo^ 

1      V  lALgl     +  vIaL^I      . 
Hence, 

AF^  <  i^-^^v)M^  +    (M2V-r'i^)!ALil 

Now      M^V-M^    =   M2(V-C)    <  0  and  1-K^V>0      , 

SO 

AF      <    (1-M,V)AL        . 
0    —  C 

Therefore,   by  Lemma  4.1 

aFq   <     -   (1-M2V)Nqj     -  {A. 2b) 

The  arguments   for  F-,{J)   are  similar.  t3 

Lemma  4.7     Under  the   conditions  of  Lemma  4.6, 

gj2)-^m(Jl)    <-Cl^j   ^Ckly..l  (4.26) 

whenever  A.,   does  not   intersect   the  boundary  that  lies  opposite  x  =  m  . 
Proof.      It   follows   from  Lemma  4.2   that 

In..  -  N..|    <  CkU.J  (^-27) 

Thus,   Lemma  4.6   implies  that 

^m(^2)   -    ^m(Jl)    <-ClNi3-Ckly.jl  (4.28) 

Note  also   that  by  Appendix  A 

2     ^  ,2 

^oj  -  "^oj     ' 

and 
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Thus,    if   ^'ijO  (x  =   0}  ^  0,    then 


F,(J2)   -    F^(Ji)   =   F^l'j^)   -   F^{Ji) 


C,N    .   +  Ck  1y    ■ 
1   oj  "oj' 


and  if   Aj -nix  =  1}     /  0,   then 


F,(J2)    -F,  (J,)iF,(J3)-F,(0,) 

Next,   let  m  =  0  and  0   <  i    <  N.      (The   argument  for  m  =   1    is  similar). 
Since  L(J)    <_  C   for  all    I-curves,   J,    it   follows   from  Lerrnia  4.2  that 

|F^(J2)    -    Fo(32)l   ICMy^jI    .  (4.29) 

The   result   (4.6)    follows   from   (4.28)    and   (4.29).  a 

Proof  of  Theorem  4.1.      Let  J(0)   =J<J<      ...<Jbea  maximal 
0  1  —      P 

sequence  of  irmediate  successors   such   that   A.  .  H  (x  =   1)    =  {5   [resp. 

t,  .  .  C\{  X  =   0}  =   0],   and  let  Y      [resp.    Y    ]   be   the  set  of  diamonds,  A  .  .    , 

I, 
crossed  by  the  above  sequence.     Note  that  t.  <   -  (1    -   x.)    if  and  only  if 

J  ~  h  1 

A  . .    c  Y      .      It   fol lows    from  Lemma  4. 7   that 
ij  0 


A..rY^        j  =  l  'if   0 


=  [F(j  )  -  F(j  )]  +  Ck    y     19.. 1 

A-  -CI 
IJ       0 

N 
By  the   results  of  Section   3,    I      \  y^- -I  £V  and  F  (J)  i  C  for  all    I-curves,  J. 

i  =  0 
Thus, 

^li;ui-K^o  • 
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(4.30) 


so  by   (4.30) 


y       N.  .   <  C 


By  a   similar  estimate  using   F-|(J)   we  obtain 


■  y      N.  .  <  C 

A  Y      ""J  ~ 


Hence,  we  have   that 


(4.31) 


(4.32) 


^  k  ^■i  <  I 

'3^'w  ''    N-j     0 


N. .  +     y        N. .  <  c    . 


(4.33) 


Mj  -■! 


Similarly,    since   V.    <  V     <   V   for  j   =    1  ,2    .    .    .  ,  we   have   that   for  every 


q  =  0,  1,  2 

such  that 


,    there   exists   a   constant   C,    indeiiendent  of  h,    k,   a,    and  q, 


I 


N.  .    <  C 


^<   t.    <  -^"tli 


i^J  - 


(4.34) 


2h        ^j  2h 

Thus,   the  estimate   of  Theorem  4.1    follows   since 


I       N..    <  qC  =   C(T) 


'3   --  ' 


ij 


for  any  q  eH   satisfying  ^  £:  T   .      This  completes   the  proof  of  Theorem  4.1 


2h 
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5.  Regularity  of  the  Approximate  Solutions 


In  this  section,  we  show  that  the  approximate  solutions,  u.  ,  are 

il 


1 


L  continuous  in  time  and  space  to  within  an  error  dominated  by  the 

mesh  length.  These  results,  stated  in  Lemmas  5.1  and  5.2,  are   needed  for 

the  proof  that  the  constructed  solution  to  (1.1)-(1.3)  actually  takes  on 

its  boundary  values  in  the  L  sense. 

Lemma  5.1  There  exists  a  constant,  C,  such  that 


Uj^(x,T2)-u^(x,T^)Idx  1  CLIt^-tJ  +  k]  .  (5.1) 


Proof:     Suppose  that 


Ol  Vl   ^^1    ^   ^m^    •••    -^   ^m'   ^2^  V-1 


(5.2) 


Then 


Iu^(x,T2)-i'^,(x,T^ )  |dx  < 


u^(x,t^-)-Uj^(x,T^)|dx 


M-1 
+  I        I    |u^(x,t.^^-)-u^(x,t  +)|dx 
J  =  m    ^ 


3=m  r 


Uy^(x,T2)-u^(x,tj^+)|dx 


u^(x,t-+)-u^(x,t.-)|dx 


0 


(5.3) 


For  s,t  c[t.,    t.      ]  we  have   that 

1 


I    |u^(x,s)-u^(x,t)Idx 
0 


ujx,s)-u,  (x,s)  Ids   + 


u^(x,s)-u,^(x,t)idx 


(5.4) 


u^(x,t)-u^(x,t)!dx 
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But 


X. 

1 


i-1 


Uu(x,s)-u, (x,t)|d) 


Thus 


£  Ch(^,2,_.|  .|v,lj) 


j  |u(x,s)-u^(x,t)Idx  <  ChL{J(j)) 


(5.5) 


Moreover,  by  (3.3)  it  follov;s  that 


max  I  u^(x,s)-u^(x,s)I  1  Cls-t^l  for  s  cJ.^^ 


(5.6) 


so  that  by  (5.4) 


u,  (x,s)-u,  (x,t)|dx  <  ChL(J(j))  +  Ck 


(5.7) 


Furthermore,  from  (5.6) 


|u^(x,t^+)-u^(x,tj-)|d> 


H  +  1 


=  J    l^h^^^^j*"'  tj-)-u,/x,tj-)|dx 


^-1 


lCh(|Y^.,^j_^  hlYi,j_il)  +  Chk   , 


1? 


so 


I 

j|u^,(x,tj+)-u^(x,tj-)|dx   <  ChL(J(j-l))   +  Ck. 


(5.8) 


Now  putting   (5.7)   and   (5.8)    into   (5.3)   and   summing  yields 


|u^(x,T2)-u^(x,T^)|dx 


<  Ch(M-rTH-l)VQ  +   Ck(M-nr+l) 
1  C[|t2-Ti    I    +   k] 

where  we  have  applied  Theorem  3.1    to  estimate  L(J(j)).     This   completes 
the  proof  of  Lemma   5.1.  3 

Next,   we   prove  a   similar  but  more   difficult  estimate   involving  L 
continuity   in   space. 

Lexnma  5.2.     Assume  that  a  =    (a,,a„,nt       ..)    is  equidistributed,   T  <  •", 
and  y-j  .y^  i  [0,1  ].     Then   there  exists  a  constant,   C  =  C(T)    such  that 


Iu^(y2.t)-u^(y^,t)ldt  <  Cly^-yJ 


(5.9) 


for  k  sufficiently  small. 

We  shall  need  the  concept  of  approximate  characteristics  to  prove 
Lemma  5.2.   It  is  convenient  for  us  to  modify  the  definitions  given  by 
Glimm  and  Lax  [2]  and  by  Liu  [6].  We  shall  define  the  set  of  approximate 


cha 


-J 


racteristics    u"    }  by   induction   on   J,   J  cH    .      An  approxi 

9    (t  f  ^i  , 


mate 


characteristic  r^   is  a  map 
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r'^  :    [t^    ,tj]   -  [0,l]x{l,2}   X  R 


where  ^^(t)   =   (x.{t),   ^.(t),y   (t))    indicates   that   the  characteristic  at 
time  t   is   at  position   x   (t),    is   travelling  on   the  £.(t)   characteristic, 
and  has   signed   strength   Y,(t).      Here   x   (t)    is  a   continuous,    piecewise  linear 
function,   and  lAt)   and  Y,{t)   are   pieccwise   constant  functions  with   respect 
to  the  partition    {t-    ,    t-   ^i,...,tj}     . 

At  fixed  J",  (f^eMj- indexes  the   characteristics   in    [0,tj],   and  Mj- 
depends  on   J"  due   to  both   the   creation  and   the  splitting  of  characteristics 
at  t    >  ty  .      To  begin   the  definition,    first   let  J  =   0.      In  this   case 
set 


M 


{0,...,N}    X   {1,2}   , 


an 


d  for    (j;  =    (i,^)cW   ,    set  Ja   =  0  and  define 


Next,   suppose  h\j_-^ 


M^  .    andfr'J'b, 


are   defined,   and  define 


J-1 


M,   =    {(}.„   ■.<t>e    M.    T,    Ge    {L,R}}U   {0,...,N}    x    {1,2}     , 
J  6  <J-  I 


with 


=  j^     for   ^,  c  Mj.T 
=  J         for  <'   =    {i,£} 


We   first  define  r'''  (t)    and  r      (t)    for   t  .  [t.    ,t   ]   and  4.  e  Mj  , 


Thus,   assume 


*R 


't'^'J-^^  ^  (>^i'^'~^)- 
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We  define  l\  (t).  Let  u^  y^   j_^    [resp.  vj^  y^  j-j  ]  denote  the  signed 

6  ' 

strength  of  the  1-wave   [resp.    2-wave]   that   lies    to   the   left  of  the 

point  (x-    1   +  "jh,tj)    [resp.{x-   +  "j^utj)]   in   the  solution  of  the  Riemann 

problem  which   is   posed  at   (x-,    t  j-, )    in   the  approximate  solution   Ul^ 

1  ? 

(so  that   (l-yi)Y  [resp. (1-H2)y^    tt]   is   the  signed  strength  of  the 

1 J  J-  I 
1-wave   [resp.    2-wave]  that  lies   to  the   right  of   (x--,   +  "jh,t-) 

[resp.    (x-   +  cijli,tj)]   in   this  Riemann   problem  solution,   0  l  Ri  ,  ^2-  ^)" 

Now  if  £=1  ,    define 

(x^(t),   l^it),   vi-,7^(t)),  t  <  tj_^,   0  <   i    <  N, 
'i^'^   =    i   (NT(^-V-l)'^'h^(^J-l)^'    t  c  Jj,   0<   i    <N, 


\ 


(0,   2,   -Vi-,Y.(tj-.i)),  t  c  Jj,    i   =  0, 


(x^(t),£^(t),    (1-Pi)7^(t)),    t    <  tj_i.   0   <   i 
r^  (t)    =    J    (x.,    1,    (l-Ui)Y.(tj_-,)),t  c   Jj,0   <   i    <  N, 


*R 


^(t-tj_^).    2,  -(l-ui)Y^(tj_i)).    t  e   Jj,    i=0    , 


and  if  £=   2,   define 


(x^(t),f_^(t),   P2Y^(t)),      t  <   tj-_^    .  0  1  i  1  N, 
r^  (t)   =     ^    (^'   ^'    VJ2Y^(tj_^)),     tcJj,   0<i<N, 


^ 


(l-^(t-tj_i),   1.   M2Y*(V-l))'    t  t   Jj.,    i   =   N   , 


(5.10a) 


(5.10b) 
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(x^(t),£^(t),    (l-v.2)Y^(t)),t  <tj_-,    ,0<i    <N, 

r^^(t)    =      <     (X.   ^^(t-tj_^),    2,    (l->'2)%(tj-l))'t^^J'   0£i    <N, 
(1,1,     (l-l'2)Yjtj_^)),     t    e:    Jj,    i    =    N    . 


The   fact  that  £    (t)   changes   value  at  t  =   t         when   i   =  0  or   i   =  N,   accounts 

for  the   reflection  of  characteristics   at  the  boundaries   x  =  0  or  1,      Morever,  when 

1-characteristics  are  reflected  at  the  boundary  x  =^  0,    the   sign  of  y   (t) 

changes   due   to   the  way  in  which  waves   reflect,   but  for(}.E;MT->  Ya,(^)   ^^^ 

a  constant  magnitude  which  we   denote  by    Iy.I    • 

Finally,  when  <^  =   (i,£),  we  have  t,   -   tj-,   and  we  define 


rJ(tj)=  (x,,£.  Nfj). 

This  completes  the  definition  of  the  approximate  characteristics.   For 

convenience,  we  let  r^  refer  to  both  the  function  defined  above,  as  well 
as  the  piecewise  linear  curve  in  xt-space  given  by  the  graph  of  x  . 

The  following  three  lemmas  give  the  important  properties  of  the 
approximate  characteristics. 
Lemma  5.3  Let  T  -  tj-  •   Then 


(5.11) 


4>  eM 


I  Iy,i  =  iui\\i\<^\i^^^), 


i,-£ 


J<J 
i,^ 


(5-12) 


and  for  j  <  J"  , 


^=^.(^3^' 


(5.13) 


the  latter  sum  being  over  all  (J"  such  that  t,  <  t-,  x.(t-)  =  x-  ,  and 

Proof.  Statement  (5.13)  follows  from  the  method  of  partitioning 
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wave   strengths    in   (5.10)   and  (5.11),   and  can  be  proven  by  a   simple   induction 

argument.      For  example,   assume  by   induction   that   the   sum  of  the  strengths 

of  the  characteristics   that  enter  a  diamond  a.,   and  have  ^,(t-)   =  t    , 

is  equal    to   the  sutii  of  the  strengths  of  the  -f-waves   th.it  enter  A..    . 

Then   the   sum  of  the  strength  of  the  characteristics   that  leave  A.-  and 

have  -^.(tj   ^  ly   must  be  equal    to   the  sum  of  thie   strengths  of  the  associated 

characteristics  that  enter  A..    ,   plus  N..,the   strength  of  the  characteristic 

'J  I  J 

that   is   created  at   (x-,t-).      By  the  definition  of  rT  . ,    this   latter  sum  equals 

I 
Y--,    verifying   (5.13).      Similarly,    (5.12)    follows   by  an  easy   induction  argument. B 

The  next   lemma    follows   immediately  by  our  choice  of woin   the   definition 

of  r       once  we  note  that,   by   (4.1), 

X   =  min   I  A  .  ( u,  )  I    >   0    . 
i,h 

Lemma  5.4     Let  iicM^  ,  j,    <  j  <  J"-l  ,  and  assume   \y Af  0.      Iflit-)   =  1 
and  a.,,    >   1    -  ^  ,   then  x'(t)   =   -  ^     for  t  e  J.^,    .      If  ^it.)   =  2  and 

"j+1    ^  ^  ,   then  x;(t)   =  ^for  tej.^^    . 

The  following  theorem  is  an  easy  consequence  of  the  definition  of 
equidistributed  sequence.  A  sequence  a=  (a-i  ,02,00, -••)  >  ctj^(0,l),  is 
equidistributed   if,    for  any  0  £  a   <  b  j£  1    , 

Urn    N([a,b],J)   ^  ^_^ 

where  N([a,b],J)  =  Card{j  <  J:  a-E[a,b]} 

Theorem  5.1  Let  a  be  any  equidistributed  sequence,  a-  c  (0,1)  . 

For  natural  numbers  m  <  n,  and  for  any  0  <  a  <  b  <  1,  let  A(m,n)  =  {m,m^-l , . . .  ,nH-nl 

and  define 
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N(n,[a,b],m)    =   Card{ j  t  A(ni,n)    :    a     E[a,b]}    , 
N(n,[a,b],J)   =:  Mjn  {N(n.  [a,b].rn) } 

Then   for  any  we  (0,1 ) , 

J"  ->  0  wJ" 

where   [x]  denotes   the  integer  part  of  x. 

Now  let  T  >   0  be  arbitrary,  and  let  J=  J(k)   satisfy  tr  -,  <   T  <  t  y  . 
Lemma  5.5     Let  ex  be  a  given  equidistributed  sequence.      Let  0  £  y-,   <  yp  i   1    . 
Then  there  exists  a  k^  such  that,   if  k  i  k   ,   then 

where  T,    denotes  the  total   amount  of  time  that  an  approximate  characteristic 
TT^with  nonzero  strength  lies   in   [y-j-h,  y2+h]x[0, tj]. 

Proof:     Let 

Myo-yJ 
^^-Tat— 
and  without  loss  of  generality,   assume  T  large  so  that  w  <  1.     Since 

tj-_i  j<  T  <   t-p  ,  we   have  J  =[t~  +   1],   so   that  we  can  make  J  arbitrarily   large 

by  choosing  k   sufficiently  small.     Thus  by  Theorem  5.1,  we  can  choose  k 

y  ~  V 
so  that  if  k  £  k   ,   then  h  <  — ^ — -     (recall   that  -r-  is  fixed),   and 

N([wJ],[0,^],J)  ^^ 

~wJ  -  ^  IT 

rind 

N([wJ],[l-  -^-  ,1],J)      ^    ^i. 
wJ  -  '^  T 

This  implies  that,  if  k  <  k  ,  then  in  any  [wJ]  consecutive  time  steps  in 
[0,tj],  we  must  have  «•  <^  -r-  l^^^P-   "-;  ±   1-  "tt]  in  cit  least  [2-1'  wJ]  of 
those  time  steps.   Thus  by  Lemma  5.4,  during  any  [wJ]  consecutive  time  steps 


.38 


in   [0,tj],   an   approximate  characteristic  with^nonzcro   strength   has   a   speed 

h  ,.„  _  ,,....  r>,k  ..r,  ^^Jy2-y^^' 


of  magnitude  ^  in   at  least    [^^  wJ]  "[-    ^  "  '~J  ^^  those  time   steps. 


Therefore,   if  k  £  k^,   a  characteristic  with  nonzero  strength  must   pass  between 


y,-h   and  y„+h  through    [y-^-h  ,y^+h']x[0,tj]   in   a  maximum  of 


,^2-^1 


xk 


[wJ]  time  steps 


But  since  an  approximate  characteristic  has  a  maximum  speed  of  y-  ,    it   can   reflect 
off  a  boundary  at  most   1  +  tj  r-  times,   so   it   can  pass   through    [yi -h  ,yp+h]x[0,t  y] 
at  most  2+tj  j^  times.      Thus   the  total    time  T     that  an  approximate  characteristic 
with  nonzero  strength  spends   in   [y,-h,  y2+h]x[0,tj]  must  satisfy 

1  C(T)(y2-yi). 

4(y        y     ) 

Note  that  k.  depends  only  on  a  and  w  = — --     .   H 

"  AT 

We  now  turn  to  the  estimate  of  Lerma  5.2.  We  first  estimate 


Uh(y2.t)-u^(ypt)|dt 


rJ 


Uh(y2't)-u^(y^,t)|dt 


and 


J-  U 


-   I 


Uh(y2>t)-u^(y^.t)|dt 


J-1 


%(y2>^'>-%(y]'^)\^^    1 


j-i 
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u^(y2>t)-u^(y-|  ,t)|dt 


^J-1 


?J 


t . 


5u^(y2,t)-Uj^(y^,t)^-[u^(y2,t)-u^(y^,t)l)|dt 


J-1 


However,   for  t  e  J-   , 


K(y2>t)-u^(yi,t)]-[G^(y2,t)-u^(ypt)il 
<   Clt-tj_^||u^(y2,t)-u^(y^,t)|    , 


and 


|u^(y2,t)-u^(ypt)i    <  C 


.1     , 


M+1         '  ^ 

hiji     +  .1  \y\j 


i^tn 


/ 


where   x     <  y,    <  x     ■,    <■  ■  ■<  y^  ±  >--   •      Hence  vie   see   that 


m 


V 


Uj^{y2,t)-u^(ypt)|dt 


(5.14) 


J-1  /    M  ,  M-1      „ 

j=OVi=nrfl  i=m     ^ 


Now  by  substituting   (5.13)    into  the   right   hand  side  of   (5.14),   we  obtain 


J-1  /  M     T      M-1   ^ 
j  =  0  \i=nr^l   ^     i-m   '-^ 


(5.15) 


I    T 


e/.l 


*  1^*'  ' 


J 


J 


where  T  is  the  amount  of  time  that  a  characteristic  r  with  nonzero  strength 
(J,  <}) 

^s  "•"  l^^n'>Vl^^f°'^J^'  ^"d  ^^1'"^  ^^^  strength 
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of  that  characteristic.      But  by  Lemma   5.5,    if  k  £  k     v/e  can  estimate 

T^  1  C(T)(y2-y^),   and  by  Lemma   5.3,       I        \y    \   £C(T),   so  we  can  estimate 

*  e  Mr 
(5.15)   by 

J-l/M  M-1  \ 

k     n  I       K^    I     +     I   Iy-J)  .<  C(T)    ly^-yj      .  (5.16) 

j=0\i=mn      ij  i=m     ^V  ^      ' 

Therefore,   substituting   (5.16)    into    (5.14)    verifies   Lemma   5.2. 
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6.      Convergence  of  Approximate  Solutions 


For  (^E  Cj^((OJ)x(0,o:,)),  we  define 


D(a,h,cj) 


[u^^.^  +   F(uj^)(^^  +  H(u^)<^]dxdt 


(6.1) 


From  (3.3)    it   follows  that 


u^(x,t)   =  U|^(x,t}f 


t 


H(u(s-t    ,    u^(x,t))dx,  te  J 


(6.2) 


Thus, 


'j-1 


u,  <;'.    dxdt 


>'   ' 


t.  0 


u,  (J,     dxdt 
h  ^t 


(6.3) 


[       H(u(5-t.,  u,  (x,t)))ds]*   dxdt 

d     f . 


Now 


u,  rf,,    dxdt 
h^i 


J 


F(u,)<>     dxdt 
n   ^x 


(6.4) 


[u^(x,t^+)<f{x,t^)   -   u^(x,t   -^)(j(x,t       )]dx 


and  from  integration   by   parts   and   (5.2) 


'o+l 


'j    ° 


[ 

J 

0         1 


H(u(s-t.,   u^(x,t)))ds](^^  dxdt 


H(u(s-tj,  u^{  x,tj^^)))ds]dx,tj^^)dx  + 
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dt 
0  t.      t. 


J   ^  [   H(u(s-t^,  u^(x,t)))ds]^  dtdx 


[u^(x,tj^^-)-u(x,tj^^-)]^(x,tj+i)dx 


(6.5) 


^J^l 


H(u,  (x,t))<j  dtdx 


"'i 


t.^,t 


f  d 


0  t 


^^  H(ii(5-t.,  u^(x,t)))<^(x,t)dsdtdx 


J   J 


Hence,  from  substituting  (5.4)  and  (6.5)  in  (6.3)  we  obtain 


Vi' 


'j  ° 


u,  <■  dxdt 
n  t 


'd.i' 


I  tf(Oh)«^  +  H{U|,)«3  dxdt 


^J   ° 


[u^(x,tj+)<>(x,tj)-u^(x,tj^^-)(;,(x,tj^^)]dx 


(6.6) 


^J^l 


-^  H(u(s-t.,  u,  (x,t)))*(x,t)dsdxdt 
dt         J    n 


^i   0   t. 


Thus,  after  summing  over  j,  v;e  obtain 


D(a,h,,f)  =  I      D.(a,h,c^)  +  I      E.(n,h,(j) 
j=1  ^         j-1  ^ 


(6.7) 


where 
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D.^   (n.h.ctO   =      [u^(x,tj^^  +  )-u^(x,tj^-,-)]4'(x,tj^l)dx 


(6.8) 


and 


Ej  +  -|(a,h,4>)    = 


J  +  1 


[F(G^)-F(Li^)]4^dxdt 


rJ  +  1 


—  H(u(s-tj,    G^(x,t)))<t  dsdxdt: 


^j     0     ^j 


(6.9) 


Now 


N        ^ri 


Dj+i(n,h,<^)    =      I        J         [u^(x-+aj^i    h,    tj^-,-)-U^(x,tj^^-)](t)(x,tj^^)d> 

^-1 


so   by   (4.16) 


D._^,(a,h,«)l    <  C(   y    \y^      |)hU|^  <  C   VhU> 


(6.10) 


Hence, 


VTh 


j  =  l    ^ 


where   the  support  of  4)   is   in   (0,"l)x(0,T) 
We  also  have   that 

1 


D.+  .(a,h,*)d^.^^|    1C(    y    I    /       1     )h2|* 
^  ^  isJi       i,j 


(6.11) 


(6.12) 


<     CV  h^U    I       =   Ch^  I  n 


X  '«> 


Hence,    for  0   <  j  -j   <  J2 


D.    {«,h,<j.)D.    (a,h,<}.)da| 

J]  J2 


(6.13) 


44 


[   D.  (a,h,(|.)da.  ]D   (a,h,<f)ciS 
■'2         2   1 


±^i'    ULI^I. 


and  where  a-  denotes  the  sequence  obtained  from  a  by  deleting  the  jp  th 
J2  -^ 

element.  By  (6.10), 


D.(a,h,(ffd:x  _<  Ch^|<,| 


(6.14) 


Thus, 


{lD.{a,h,<t>)fda 


(6.15) 


21  . 


D.  D.  da  +  5; 

J]  J2       j  J 


2 
D.da 


lilfc^yULKl, 


+  ({)^cv  ^h^\<^\^  =  ch(U|^  +  U  1^) 


X'co 


So,    for  each  4.cC   ( (0,1  )x(0,"') )    there  exists   a  sequence  h      »■  0  and  a   set 
Ac.g   [0,1]  with  measure  one  such  that  if  acA,    then 


D.(a,h    ,4.)  ->  0     as  h      >  0 
J         m^'  m 


(6.16) 


Next,   let   (O   )       ,,  be  a   countable   dense  set   in   c"(  (0,1  )x(0,-) )   with   respect 
n   nc  M  0 

to   the  norm   |      |^   .     Then  after  refining  h     above  by  a   diagonal   argunsnt  and 
after  taking  the  appropriate  countable   intersection  of  sets  A     of  measure  one 
we   can  conclude   that   there  exists   a   sequence  h     ♦  0  and  a  set  A^.5i[0,l]  with 
measure  one   such   that   ifacA,   ntN^    then 
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I   ^i(ci,h|^,<^^)  -.  0  as    m  -^oo  . 


('"^.17) 


It  follows   by   (6.11)    that    acA,  (t.  e   C^   ( (0,1  )x(0,") )  ,   then 

I  D^(a,h^,*)   -  0     . 


(6.18) 


Now  by   (4.10) 


t.   .    1 

J+1 


[F(u^)   -   F(u^)]6dxdt 


t.        0 
J 


<     Ck' 


(6.19) 


Also,   for  X  el  •  , 


^3-1    ' 


^H(Li(s-t^,  L^(x,t)))^.dsdt 


t- 


Vi  b^i 


t.       s 

J 


<     Ckl 


H     u     u,  ,   <.dtds' 
u     u     nt 


1    ivi" 

t. 
J 


(6.20) 


1     CkUl^dy  I    + 

i-l,j 


l^lol' 


Hence, 


tj.l   '     ' 


.      0    t. 


^H(ij(s-t^    ,   u^(x,t)))*dsdxdt 


<  C(   I    IyI  J)khl4> 


(6.21) 


1  C  V  khUL  1  CkhUI 
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So,    it   follows   from   (6.19)   and   (6.21)    that 


I    I   E    («,h,^)  I  <   CTkkl     +   CThUI 
j=l    J 


(6.22) 


Thus,    from  (6. 1  7)  and  (6.22)we  can   conclude   that  if  acA  and   (}icC'°(  (0,1  )x(0,oo) ) , 
then 


D(a,h    ,  <))   ->  0  as  171  -»■  «>. 


(6.23) 


Next,   choose  T  >  0  . 

Let    {tg}  C  (0,T)   be  countable  and  dense  and  letacA  be  fixed.      Since 

Var       u   (•  ,t)   £  C 
xr[0,l] 

for  te(0,«=),   Helly's   theorem  implies   that   there  exists  a   function 

u(-,t):    {t„} -^L    (0,1)    such   that   (after  further  refining    {h^})    for  pcB 

Ul,    (•,   t   )   ->■  u(-,t   )    in   L    (0,1)   as  m  ->-  CO  .      it   follows   from  Lemma   5.1    that 
"m  P  f^ 

u,     (-.t)   converges   to  u(-,t)    in  L    (0,1)    for  each   tt:(0,T)   and 


|u(.,T   )   -   u(-,T,)|      ,  <  C|t^   -    T, 

'^  '      l'(o,i)  ~       ^        ' 


for  T| ,    T2  c(0,T).      Hence: 


^%   *t  ""   ''^%   ^*x  ^   "(^h   )*3dxdt 


0  0 


m 


m 


m 


(6.25) 


[u*^  +   F(u)4.^  +   H(u)4']dxdt 


0  0 


as  m  -*■  0 


Since   D(a,h_,({))  ->■  0,   u  is  a  weak   solution  to   (1.1).      Furthermore, 


m 


since   u^(0)   ->   u^   in   L    (0,1), (5. 24)    implies    that 
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Urn     u(-  ,t)   =   u^  (6.26) 

t  ^  0 

To  obtain    (1.13)   note   that  the  equidistributed  sequences  have  measure  one 
in    -n, [0,1],    so  that  we  can   further  assume  that  «   is  equidistributed.      Since 

u^    ^    u   in   l\(0,1)x(0,T))    , 
m 


we  also  have 


u^   (x,-)  -^  u(x,-)    in   l\o,T)  (6.27) 

m 


for  almost   all    x.     An  application  of  Lemma  5.2  now  implies 

u,    (x,- )  -+  u(x,- )    in   I.   (0,T) 
m 

for  all    x  and 


|u(y,,-)    -   u(y   ,.)!    1    ,  <   C|y     -    y^l 


But 


(6.28) 


P.     (0,-)   ->  pg   in  l\o,T) 
"m 

so   that  we  can   finally  conclude   from  (6.28)    that 

^OTi     p(x,-)    -  p„    ,  Urn       G(x,-)    =   0   in   l\o,T).  (6.29) 

X  -  0  '^       X  -V  1 

This  complete  the  proof  of  Theorem  1. 
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Appendix  A 

In  this  appendix,   we   give  the   proof  of  Lenirna  3.3.     He   relax   the 
assumption   that  H   is   locally  Lipschitz   at   [Gj    =   G     so   that  v^e   can  better 
model    the   transition   from  laminar  flow  to   turbulent  flow.     We   replace 
property   (1 .8)   by: 

H   is  Lipschitz   continuous   on    [-G   ,   G  ]  and  H   is  locally  Lipschitz         (A.l) 

continuous   on    (-^-.-G   ]u[G   ,«>). 

Since   (1.7)    is   to  remain  valid  at    |G|      =   G^   in   the  sense  of  distributions, 

we  have  that 

Hi-G  -)=  Urn       H(G)   <     tim      H(G)=H(-G+) 
^       G  +  -G^  G4-G^  ^ 

(A. 2) 

and 


^^(C,^-)-  tljv        H(G)  ±lXin        H(G)    =    H(G^+) 


G  +  G^  G  4  G^ 


For  the   remainder  of  this  appendix,  we   shall    consider   H  to  be  multi-valued 


at  +  Gj,  so  that 


and 


H(-G^)   =    [//(-G^-),    H(-G^+)] 


For  G  /  j^  G^,    H  remains   single   valued.      Note   that  since  H    is  a  maximal 
monotone  function,    the   differential    equation 

(where  H(u)   =    (0,   -H(G)/p)        is  now  a   multi-valued  map)   always  has  a 
unique   solution    [7].      Therefore,    the   friction   step   (3.2),    (3.3)    is   still 
wcl 1 -defined. 

Now  let   F(r,s)   h   -H(G)/p.      Only   the   following  monotonicity  property  of   F 
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shall    be   used   in   the   proof  of  Lemma   3.3: 

^  {r^s)  ^    F(r2,s)  when   r-^  >    r^,      Ivl<    c    ,  (A. 4) 

F(r,s^)l    ^(r.s^)   when   s^j  >    s^,      |v|  <    c    .  (A. 5) 

When    F  is  multi-valued  at    (r,s),   we   take   this   to  mean   that   (A. 4)    and   (A. 5) 
are  valid  for  every  element  of   F(r,s).      The  monotoniuty  assumptious    (A. 4) 
and   (A. 5)   are   consequences  of   (1.7)    since 


dr  -  -  p    lngvi-^v   -  G     c- 


-!7'=-^[H^(l4)-^^]lO      , 


for  I  vl  <c  .   Given  these  assumptions,  we  now  state  and  prove: 
Lemma  3.  3  If  J^  is  an  immediate  successor  of  J-,,  then 

Lg(J2)  +  L(J2)  1  Lg(J2)  "^  "-(J^)   .  (A. 6) 

Proof:   First  consider  the  boundary  cases  i  =  0  and  i  =  N  .   It  is  easy 

to  see  that 

2    ~2 
oj     oj 


and 


so  (A. 6)  follows  for  these  cases.  Now  let  L[U|  ,Un]  denote  the  sum  of  the 

strengths  of  the  waves  in  the  Riemann  problem  solution  [u,  ,Un].  Also,  let 

u,  (t)  =  u(t,u,)  and  Un(t)  =  u(t,Un)  .  To  verify  (A. 6)  it  suffices  to  show  that 

LCi^L^^)'  ^r(^^^  -  L[\'^'r^'         t>0.  (A. 7) 

So,  consider  the  case  sign  v.  f   sign  Vp  .   Note  that  since  h(0)  =  0,  we 

have  that  sign  Vj.(,;)  /  sign  v^(c;)  for  some  ^  e  [0,°^)  if  and  only  if  signVp(t)  / 
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sign   V[^(t)    for  all    t  e    [0,^).      In    this   case,    (A. 7)    follows   since 

In  P|(t)   =  In  pp(t)   and   since    |vp(t)    -    V|^(t)|    is   decreasing   for  tf  [O,'-). 

For  the  other  case,   assume   sign   V|^(t)    -    sign   Vn(t)    >  0   for  all    t  c    [0,°°) 
(the  case  v,_(t),   Vf^(t)    <  0   is   similar).     Refer  to  Fig.    A.l    to  sec  that 

if  u   (t)    lies    in  quadrant   II   or  quadrant   IV  with  respect  to  u   (t),    then 
L[uL(t),    Up(t)]  -    \Jin  PL(t)    -  Zn  p^{t)\ 
=    \ln  p^   -  In   p^\    <_  L[uj_,Uj^]      . 


cZnP  1         V 


>r 


II 


III 


u^Ct) 


IV 


Fig.    A.l 

Finally,   suppose   that  Up(t)    lies    in   quadrant   I    [resp.    quadrant      III]  with 
respect  to  U|^(t).      By   (A. 4)   and   (A. 5)    it   follows   that  if  Uj    [resp.    Ujjj] 
is   in  quadrant     I  [resp.   Ill^with   respect   to   u,  ,    then 


|H(u„i)l   1   1H(ul)1    i  |H(ui) 


(A. 8) 


Thus,  we  conclude  that  Un(c)  is  in  quadrant  I  [resp.  Ill  ]  for  all  earlier 
times  CE[0,t].   But  by  (A. 8)  Vp(c)  -  V|_(rJ  >  0  [resp.  w^{r,)   -    v^(0  >  0] 
for  all  times  c  e   [0,t].   This  implies  that  L[U|^(r),  Up(:,)]  is  decreasing 

throughout  [0,t],  so  that  (A. 7)  follows  for  this  final  case.  □ 
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Appendix  B 

In    this   appendix,   vie   give   the   proof  of  Lemma  4.4. 
Lemma  4.4      If  0   <   i    <  N,    then 


aR"^!    <  }  D^j,  -e  =   1,2    .  (B.l) 


Proof.  We  do  the  case  £  =  2.   It  suffices  to  consider  the  case  when 

only  two  1-shocks,  S:  and  si    ,   enter  A..  [13].   In  this  case,  a  1-shock, 

1  2 

S       ,   and  2-rarefaction,   R     ,    ,    leave   a..    .     Without  loss  of  generality, 
out  out  ij 

assume    |SJ    >_  \S    \    .      In  order  to  obtain    (B.l),   wo   first   must   estimate 
derivatives  on   1-shock   curves.      So,    if  v  -   v(z)   is  the   velocity  along  the 
1-shock   curve   that  starts   at   u    ,  then   by   (2.9) 

,    -z/2       ^  z/2s 
v  -  v.    =■   c(e     '      -  e         )    , 

z  =  £n  p,    -  -Cn  pp   <  0   . 


Thus, 


L        -"   ^R  - 


^=-ccoshf  (B.2) 

dz  2    • 

We  estimate  4"^- with  the   following  technical    lemma: 
dz 

Lemma  B.l      If  w  ^0,    then 

■^  w  ^  l-sech  w.  (B-  3) 

Proof.     Let    p(w)    =  |  -    {l-sech  w}    .      If    p'(w)   is  nonnegative,   then    p 
is   increasing;   and  since    P(0)   =  0,    this  would   imply  that    P(w)  ^  0  for 
w  ^  0.      Now 

P'  (v;)    =  -  -sech  w  tanh  w   . 

But  it  is  easily  verified  that  the  maximum  of  sech  w  tanh  w  is  ^  > 

so  that  p'  (w)  >  0.   0 
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cX)iP 


Fig.    B.l 


Now  consider  the   interaction  of  tv/o  1-shocks   diagrammed  in   Fig.    B.l 
Here 

,2 


A  =   c   R 


out 


C  -   c  sj     <   D  =    |v(S^^^.)    -   v(s])|      , 


B  =   cIS^    t   -   s! I    =   C-A 
I    out  1  I 


{B.4) 


(recall    that  we   use  the  same  notation   to  denote  both  a  viave  and   its   strength) 
By   (B.2), 


D  =   ^il^     dv    ^^  ^pl    "   ^^^^   ^1 
where    |-j^  (2   S,  )  |    denotes    the   derivative  -.-  of   the   1-shock  curve   starting 


(B.5) 
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at  u.  ,   evaluated  at  z  =  2  S,    .      But 


A  _   T        B       ,        B 
C-   ^    -   C-^    -   D 


I   -  —  <1    -— -   1   -n 


Thus  bj'   (B.4)    and   (B.5) 


R^   4-    <  (l-n)si    <   {1-sech  s]  }  s],        .  (B.6) 

out   —  2   —  12 


By  Lemma  B.  1  , 


{1-sech   S    }    <  ^J_     >  (B-7) 

2 


so   that  substituting   (B.7)    into   (B.5)  yields 


R^    <-ls;    S?   =  iD^ 
out  ^  2     1      1        2     ij 


This  completes  the  proof  of  Lemma  4.4.  E3 
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contractor. 
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